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Abstract 

We study the dynamics of Josephson junction arrays with posi- 
tional disorder and driven by an external current. We consider weak 
magnetic fields, corresponding to a frustration / = n + 1/25 with n 
integer. We find that above the critical current i c there is a plastic 
flow of vortices, where most of the vortices are pinned and only a few 
vortices flow through channels. This dynamical regime is character- 
ized by strong fluctuations of the total vorticity. The number of the 
flow channels grow with increasing bias current. At larger currents 
there is a dynamical regime characterized by the homogeneous mo- 
tion of all the vortices, i.e. a flux flow regime. We find a dynamical 
phase transition between the plastic flow and the flux flow regimes 
when analyzing voltage-voltage correlation functions. 
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1 Introduction 



The problem of transport in a medium with quenched disorder leads to very 
rich physics in a great variety of systems Q . Examples include charge density 
waves, Wigner solids, magnetic bubble arrays, fluid invasion in porous media 
and pinning phenomena in type II superconductors p], |3|, §], |5j, ||, [7j. All 
these systems have in common a competition between interaction forces and 
randomness under an external driving force. For very large driving forces 
the randomness is believed to be irrelevant, and the dynamics is interac- 
tion dominated. The system shows homogeneous flow, and all the internal 
degrees of freedom move as a whole. On the other hand, for low driving 
forces the dynamics is dominated by the disorder. The system shows a very 
inhomogeneous flow, where the motion breaks up into pieces moving with 
different velocities or not moving at all. It has been recently discussed that a 
dynamical phase transition may exist between these two kinds of flow |3], [7| 

In the last years, a great deal of attention has been devoted to the study 
of vortex pinning and dynamics in type II superconductors as an example 
of this kind of problem [J5], || |7|. It has been shown that for all but very 
weak pinning, dissipation starts at the critical current via channels of flow 
of vortices. Experiments by Bhattacharya and Higgins |J have studied the 
different dynamical regimes of plastic flow and flux flow motion in super- 
conductors. Koshelev and Vinokur []7| have discussed a dynamical melting 
transition from a moving vortex lattice at large currents to a fluid flow of 
vortices at lower currents. 

Most of the experimental systems mentioned above have the difficulty that 
there is no control of the nature and amount of disorder. Therefore the study 
of disordered Josephson junction arrays ( JJA) ||, ^ |10[ becomes particularly 
promising here, since they can be specifically fabricated with controlled ran- 
domness Also, the recent developement of imaging techniques in JJA 
can allow for a direct observation of the different spatio-temporal dynamical 
regimes expected in these systems [ |TT| , [T2J]. In this report we consider JJA 
with positional disorder || under an external magnetic field. In this case the 
effective amount of radomness can be changed by increasing the magnetic 
flux per plaquette by integer multiples of the quantum of flux. We have 
previously shown that in the limit of very large disorder (gauge glass) this 
system has a dynamical phase transition at the critical current ||13|| . Also, 
the dynamical phase diagrams for different amounts of frustration and bias 
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currents were discussed in [[L4J]. In this contribution, we consider the case 
when the average number of quantum of flux per plaquette is / = n + 1/25, 
with n integer. In the absence of disorder, this magnetic field corresponds 
to a ground state consisting in a periodic vortex lattice. Here we will study 
the transition between the regime of plastic or inhomogeneous flow and the 
regime of flux flow. 



2 Dynamical equations 

The dynamical equations of current driven JJA can be obtained from con- 
sidering the resistively shunted junction model for each Josepshson junction 
plus current conservation at each node fllS] , p~3f| . One obtains the following 
set of dynamical equations for the superconducting phases 9(r): 

^£G(rV){/^(r) (1) 

£ h sm(0(f + A) - 0(f) - 2tt/ a (0)} , 
±fi 

where 9(f) is the phase of the superconducting wave function at site f , fi = 
%iV\ /a(^) = <jr A.dl is the line integral of the vector potential with $ 
the flux quanta. The array has N x N sites, and a is the lattice constant. 
Jo is the critical current of each junction, assumed to be independent of 
disorder. All the effect of the randomness is taken in f/x(f). In the Landau 
gauge f(f,f) = —7^ — rx ^ Ty+rv \ with H the applied magnetic field. We 
take random displacements of the sites, f/a = (n x + 5 X , n y + 5 y ), with n x , n y 
integers, and 6 x,y a random uniform number in [— A/2, A/2]. Noting that 
(X)g//i(f)) = Ha 2 /Q = f, we consider here f = n + 1/25. The effective 
amount of disorder is W = nA |J, so one can increase the disorder by 
increasing n (i.e. the magnetic field). In the absence of randomness (A = 0), 
/ = n + 1/25 is equivalent to / = 1/25. TZ is the shunted resistance of each 
junction and G(f,f) is the two dimensional lattice Green's function. The 
boundary conditions are periodic along the x-direction. At the bottom (top) 
of the array the external current is injected (taken out) with I ext (r y = 0) = I, 
(I ext (r y = Na) = —J), and I ext (f) = otherwise. We evaluate Eq. ([[]) with 
the same fast algorithm as in Ref. fl5| . The time integration is done with a 



d9(f) 
dt 
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fourth order Runge-Kutta method, with integration step At 
(tj 



0.01 -O.lrj 



2 J 1 ^/ ) ; for time intervals of t = 5000rj, after a transient of 2000tj. 



Magnitudes of interest are the voltage drops along the direction of the 
current, 

f d6{r + y) _ d6{f)\ 
dTJ 



v y (r,t) = Tj 



dt 



(2) 



and the distribution of vortices 



n(R, t) = — nint 



6{r + p) - 6{r) 
2^ 



hif) 



(3) 



with nint (a?) the nearest integer to x, and the sum runs around the pla- 
quette R. We will study time averaged quantities such as the voltage v = 
i) Yjr v y{^i t)> with v the average over time. The total vorticity n<r(i) = 

J2ii n (R, t) satisfies n T (t) = f, because of vortex number conservation, 



i 



N X (Ny-l) 



(iV-1) 2 



but it can have temporal fluctuations (Snr) = ^t(^) 2 — nr(f) 



3 Transition to the flux flow regime 

If the frustration is / = n + p/q and <C 1 there is a dilute vortex 
configuration in the array. Here we consider p/q = 1/25 in arrays of size 
50 x 50. In the absence of disorder (A = 0), the vortices form a periodic 
lattice in the ground state. For intermediate values of disorder W ~ 0.2 — 0.5, 
the equilibrium configuration of vortices at zero bias J = is a random 
array of vortices, with no crystalline order. (We find that a distorted vortex 
lattice exists only for very weak disorder W < 0.1). Here we consider the 
particular case of A = 0.05 and n = 4 (W = 0.2). A full diagram of the 
different dynamical regimes as a function of W and bias current / was given 
in Ref. JT4] . When solving numerically Eq. ([I]), we find that well below a 
critical current I c all the vortices are pinned by the disorder potential and 
there is no dissipation. Close but below I c some vortices are depinned and 
move during a transient time until they are pinned again in a new stable 
configuration. At I c there is one channel of vortices that move without being 
pinned again. Close and above J c only few vortices move (the others remain 
pinned) following certain channels of flow. An example of this is shown in 
Fig. 1(a). As the current is increased further, more vortices are moving and 
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therefore more channels of flow are opened [Fig. 1(b)]. When increasing 
the current even more, eventually all the vortices can move but still with 
a random, fluid-like, inhomogeneous motion [Fig. 1(c)]. For this range of 
currents the volt age- current curve (see Fig. 2) shows a nonlinear increase of 
the voltage v from the critical current i c (normalized by the Josephson current 
ic — Ic/Iq)- This regime of inhomogeneous vortex motion is usually called 
plastic flow, and it has been studied in disordered type II superconductors 
[|]. In disordered JJA, it can be characterized by fluctuations of the total 



vorticity 8nr > |TJ|, |TJ|, i.e. by the existence of flux noise. At high 
currents, above a characteristic current i p , above which 5nx = 0, the voltage 
v grows linearly with the current with an slope proportional to the vortex 
concentration p/q. This is the so called flux flow regime. In this case all the 
vortices move with almost the same velocity. 

A quantitative understanding of the transition from the plastic flow to the 
flux flow regime can be obtained by analyzing the correlations of the average 
voltages. Let us consider the voltages along the direction of the current, 
v y (f, t), which are a measure of the vortex velocities. The homogeinity of the 
vortex motion can be studied from the time averaged voltages v y (r) = v y (r, t). 
We define the correlation function 

N x (Ny - 1) £^(rK(r + xe x ) - v 2 , (4) 



and similarly along the y direction, T vy (y). In Fig. 3 we show the normalized 
correlation C x (x) = T vx (x)/T vx (0). We obtain that the voltage correlation 
can be approximated by the form 

T vx (x) ~ T ra (0) exp(-x/£ ra ). (5) 

Therefore, the voltage correlation length £ vx is a measure of how inhomege- 
neous is the flow along the x direction. In other words, it is a measure of 
the typical distance that takes a vortex to turn from a straight path perpen- 
dicular to the current, so that £ vx = oo (i.e. C x (x) = 1) means a perfectly 
straight motion. We see in Fig. 3 that there is a dynamical phase transition 
at the characteristic current i p (in this particular case i p = 0.405) where 
£„ x diverges. We obtain that for i > i p , in the flux flow regime, C x (x) = 1 
(Cm = i- e - an the vortices move in a straight path, with the same velocity 
along the x direction [this is also confirmed by analyzing the voltages v x (r)]. 
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On the other hand, the voltage-voltage correlations along the y direction, 
T vy (y), never reach complete homogeinity, therefore above i p the voltages 
Vy{f) have a dependence on y only. This dependence consists on small fluc- 
tuations of v y (f) = v y (y) along the y direction around the mean value v. 
This means that in the flux flow regime, i > i p , the effect of the disorder 
potential becomes irrelevant only along the direction of the Lorentz force 
(the x direction), but along the y direction the disorder is still important. 
Another interesting result is that for low currents just above i c , the length 
£ vx is almost constant, as it can be seen in Fig. 3. This corresponds with the 
regime in currents where some of the vortices flow in channels while some 
other vortices remain pinned. At high currents, of the order of the Josephson 
current Iq, vortex-antivortices are induced in the JJA. A transition from the 
flux flow regime to dynamical regimes dominated by the vortex-antivortex ex- 
citations appear at these very high currents (involving again plastic flow and 
homogeneos flow, but of vortex-antivortex pairs). These dynamical regimes 
were discussed in ||14|. 



4 Discussion 

In conclusion, we have shown that the transition from the plastic flow regime 
to a flux flow regime seems to occur as a dynamical critical phenomenon 
with a diverging correlation length. These could be studied experimentally 
in JJA with controlled disorder at low temperatures using recently developed 
vortex imaging techniques [[□]] . Another possible test of this transition from 



vortex random motion to vortex straight motion could be obtained from 



measurements of the fluctuations in the Hall voltage |16 . 
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Figure Captions 



Figure 1: Diagrams of vortex flow in 50 x 51 Josephson arrays with frus- 
tration / = 4 + 1/25 and disorder A = 0.05. (a) i = 0.205 (i c = 0.203), (b) 
% = 0.21, (c) % = 0.35, (d) % = 0.6 (i p = 0.405). The black squares represent 
the position of the vortices at the present time, and the gray squares the 
positions where the vortices have been previously. Therefore the gray lines 
indicate the paths of the vortices. 

Figure 2: Voltage-current curve for a 50 x 51 Josephson junction array 
with frustration / = 4+1/25 and disorder A = 0.05. Currents are normalized 
by I and voltages by TZIq. 

Figure 3: Left: plot of the voltage- voltage correlation function C x (x) for 
different currents. Right: plot of the correlation length £ vx extracted from 
C x (x) as a function of current. 
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